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Recent experimental and theoretical research in the area of neutrino 
interactions in the ~ 1 GeV region are reviewed including topics like: the 
problem of value of quasielastic axial mass, neutral current ir° production, 
coherent pion production. Many comments are devoted to status and cur- 
rent development of Monte Carlo events generators. 
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1. Introduction 

Knowledge of neutrino interactions in the ~ 1 GeV energy region is 
important because this energy domain is typical for majority of neutrino 
oscillations experiments performed during recent ~ 5 years and also those 
scheduled for the near future. The list includes K2K, MiniBooNE, Sci- 
BooNE, MINOS, T2K and NOi/A. The only exception is OPERA with 
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(^ higher energy neutrino beam from CERN. 



Neutrino oscillations are the energy dependent phenomenon and the most 
straigthforward analysis of experimental data requires reconstruction of neu- 
trino energy. The neutrino flux spectrum is typically rather wide-band (de- 
spite significant improvements introduced with the idea of off-axis beams) 
and the interacting neutrino energy must be estimated based on the ob- 
servation of the leptonic and/or hadronic final states. The precision of the 
analysis depends on the knowledge of neutrino interaction cross-sections, 
both inclusive and exclusive in several most important channels. If the 
analysis of the oscillation signal does not include as the intermediate step 
the neutrino energy reconstruction [T], it is still based on the comparison 
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with predictions from the Monte Carlo (MC) events generators and rely on 
how well they are known and implemented in numerical codes. In the anal- 
ysis of the oscillation appearance signal as seen in Cherenkov detectors it is 
very important to evaluate the background from NC-zr production events. 
Since the required precision of new oscillation experiments is higher with re- 
spect to what was sufficient until recently the series of dedicated workshops 
was initiated with the aim to organize a forum of discussion for experimen- 
talists and theorists, including specialists in nuclear physics, to present new 
measurements and proposed improvements in models [2J. Thanks to Nulnt 
workshops during the last 8 years the knowledge about neutrino interactions 
improved significantly even if some problems and limitations turned out to 
be quite robust. 

In the ~ 1 GeV energy region one distinguishes several processes that invoke 
quite different theoretical descriptions. The terminology is sometimes con- 
fusing because neutrino reactions can take place both on free nucleons and 
(this is the common situation) nucleus targets. Fortunately in the ~ 1 GeV 
energy region typical values of momentum transfer are large enough and 
with a good approximation one can assume that neutrino-nucleus reaction 
occurs on individual bound nucleons (impulse approximation - IA). Limita- 
tions of this picture will be addressed several times in what follows. Thus 
we distinguish: quasielastic (QE) CC reaction 

v + n — > l~ + p 

v+p -+l + + n (1) 

and its NC elastic counterpart 



v + N ->• v + N 

v + N ^v + N (2) 

single-pion production (SPP) reactions 

v+p — >l~+p + TT + 
v + n — > l~ + n + 7r + 
u + n -> l~ +p+TT° (3) 

(with many other channels for anti-neutrinos and for NC reactions) and 
more inelastic reactions (the misleading name DIS is commonly used) . The 
seperation of three types of reaction is quite general and the specific feature 
of the ~ 1 GeV region is that all three (QE, SPP, DIS) contributions to 
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the total cross-section are important. The necessity to have a consistent 
inclusion of both SPP and DIS dynamics is a source of many difficulties 
in the construction of MC generators of neutrino interactions. In the case 
of neutrino-nucleus reaction there is a possibility also of a coherent pion 
production (COH) 

v + ^X -^v + Tr° + iX. (4) 

The second reaction contributes to the overall NC ir° production and is a 
subject of intensive experimental and theoretical investigation. 

Nuclear effects make this picture more complicated. The first obvious and 
nontrivial difficulty is the problem of nuclear enviroment modifications of 
free neutrino- nucleon interactions. The second problem comes from final 
state interactions (FSI). Even if we accept as reasonable the assumption 
that the reaction takes place on individual nucleons, the particles that arise 
in the final state propagate through nucleus before they can be experi- 
mentally detected. Thus experimentalists would rather like to speak about 
QE-like events defined by the condition that there are no mesons in the final 
state or SPP-like events with only one pion and no other mesons in the final 
state. It is clear that there is a significant difference between QE and QE- 
like events because the latter include those in which pions were produced 
in the initial interaction but were later on absorbed inside the nucleus. It 
is also possible that rescattering of a nucleon from QE primary interaction 
results in pions in the final state. 

In direct experimental analysis further complication comes from the detec- 
tion thresholds for kinetic energies of various species of particles. These can 
be very different. In the water Cherenkov detectors like SuperKamiokande 
they are ~ 75 MeV for n and ~ 485 MeV for protons while in the ND280 
near detector in the T2K experiment they may be about 100 MeV and 
150 MeV correspondingly. Thus the measurement of QE-like and SPP-like 
cross-sections require also some information from MC events generators and 
the uncertainty introduced by them vary from experiment to experiment. 
FSI effects are an important source of uncertainty of MC events generators 
and thus in the data analysis as well. This is why recently the experimental 
groups tend to provide the neutrino-nucleus cross-section data (with the 
FSI effects included) rather than neutrino-nucleon ones. Such data include 
corrections for detector efficiency and the effort is done to make them inde- 
pendent (as much as possible) from the nuclear physics assumptions of the 
MC codes used in the analysis. 
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Many nuclear physics complications are absent in experiments performed 
on the hydrogen or deuterium target. This is the reason why there is still 
an interest in old bubble chamber experiments like ANL or BNL and surpris- 
ingly they still can be a source of interesting information [3j. The advantage 
of new experiments like MiniBooNE and SciBooNE is in much higher statis- 
tics and better understanding of the neutrino flux. 

In my presentation I will include many comments on the status of Monte 
Carlo events generators. It is important that the generators are modernized 
in parallel with better understanding of ~ 1 GeV neutrino interactions. 



2. Quasielastic axial mass 

Charge current quasielastic reaction (CCQE) is the dominant process in 
the sub-GeV neutrino energy region typical for MiniBooNE, SciBooNE or 
T2K experiments. The theoretical description of neutrino-nucleon reaction 
is based on the conserved vector current (CVC) and the partially conserved 
axial current (PCAC) hypotheses. As a result of simple analysis the unique 
unknown quantity is the axial form- factor Ga(Q 2 ) for which one typically 
asssumes the dipole form with only one unknown parameter called axial 
mass Ma- If the deviations from the dipole form of Ga are of similar size to 
those in the case of electromagnetic form-factors it would be very difficult 
to detect them and the basic assumptions described above seem to be well 
justified. Thus the aim of CCQE experiments is to measure the value of 
Ma- Even if in the experiments neutrinos interact with bound nucleons the 
reported results should always refer to the parameter in the formula for free 
nucleon scattering. Obviously, any such measurement done on a nucleus 
target contains a bias from the model of nucleus used in the data analysis. 

The measurements of Ma typically focus on the analysis of the shape of 
the differential cross-section in Q 2 that turns out to be sensitive enough 
for quite precise evaluation of Ma- The investigation of only the shape of 
the distribution of events in Q 2 has an advantage that it does not rely on 
the knowledge of the overall neutrino flux that usually carrries much uncer- 
tainty. The dependence of the total cross-section on Ma can also be used 
as a tool to fix its value. The limiting value of the CCQE cross-section as 
E v — >■ oo can be calculated in the analytical way assuming dipole vector 
and axial form- factors [4J. In the exact formula the dependence on Ma is 
strictly speaking quadratic but in the physically relevant region it is with a 
good approximation linear. It is an important fact that if the value of Ma 
is increased e.g. from 1.03 to 1.23 GeV the cross-section and the expected 
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number of CCQE events is raised by ~ 20%. 

In the past there were several measurements of Ma mostly on the deu- 
terium target and until few years ago it seemed that the results converge 
to a value of the order of 1.03 GeV. There is an additional argument in 
favor of a similar value of Ma coming from the weak pion-production at 
low Q 2 . PCAC based evaluation gives the value of 1.077 ± 0.039 GeV 0. 
On the contrary all (with the exception of the NOMAD experiment) more 
recent high statistics measurements of Ma report much larger values: K2K 
(oxygen, Q 2 > 0.2 GeV 2 ) -> 1.2 ± 0.12 ©; K2K (carbon, Q 2 > 0.2 GeV 2 ) 
->• 1.14±0.11[7J; MINOS (iron, Q 2 > GeV 2 ) -► 1.19±0.17; MINOS (iron, 
Q 2 > 0.3 GeV 2 ) -)• 1.26 ± 0.170; MiniBooNE (carbon, Q 2 > GeV 2 ) 
->• 1.35 ±0.17; MiniBooNE (carbon, Q 2 > 0.25 GeV 2 ) -»■ 1.27 ±0.14 [9] (for 
completness: NOMAD (carbon, Q 2 > GeV 2 ) ->• 1.07 ± 0.07 [10]). 



There are a few possible explanations of the discrepancy. In the simplest, 
one notices large uncertainties of the measured value of Ma and treats the 
discrepancy as merely statistical fluctuations (after all the effect is on the 
< 2o" level). The problem is that there are several independent measure- 
ments. In the case of MiniBooNE large values of Ma were obtained from 
the investigation of the shape of the distribution of events in Q 2 and also 
as a fit to the normalized cross-section and both evaluations do agree. This 
weaken doubts that are sometimes raised concerning the MiniBooNE's un- 
derstanding of the overall normalization (integrated flux). In fact, there are 
other MiniBooNE measurements e.g. (CCl7r + ) giving rise to larger than 
expected cross-section. The normalization (flux and fully correlated sys- 
tematic errors) uncertainty is evaluated by the MiniBooNE collaboration to 
be 10.7%. A delicate element of the MiniBooNE's data analysis is a sub- 
traction from the sample of QE-like events of those that are believed to be 
not QE in the primary interaction. In the analysis the NUANCE [llj MC 
event generator based on the Fermi gas model was used. Obviously such 
subtraction depends on assumptions of the MC model. MiniBooNE collabo- 
ration corrected the MC prediction for this background by the function that 
was obtained by comparing a sample of SPP-like events to the predictions 
of the same MC generators. The shape of the correction function is rather 
poorly understood but it has an obvious and quite important impact on 
the extracted value of Ma- The function quantifies a lack of precision in 
describing processes like pion absorption and this can have different effect 
on understanding of QE-like and SPP-like samples of events. 

A separate difficulty is related to the low Q 2 region. It has been known 
for many years that MC events generators have problems with correct re- 
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production of the shape of differential cross-section in this region (low Q 2 
means typically Q 2 <^ 0.1 GeV 2 ). This is the reason why in the data 
analysis very often (see above) some cuts are imposed. The low Q 2 problem 
has to do with the validity of the impulse approximation. We know from 
the electron scattering that for momentum transfers q < 350 — 400 MeV/c 
the models based on IA fail to agree with the data. In this region collective 
nuclear effects become important and computational techniques like RPA 
or CRPA should be used [12J. Since q > to where uj is the energy transfer, 
and Q 2 = q 2 — uj 2 , the region of the failure of IA is clearly contained in 
the domain Q 2 < 0.1 GeV 2 . How large is this dangerous region? Contrary 
to what might be expected, in the case of neutrino interactions it is always 
large, of the order of 15% — 20% of the total cross-section, independently on 
the neutrino energy (for energies E v < 500 MeV it is even larger) [13]. Ex- 
perimental groups invented some ad hoc solutions to deal with the low Q 2 
problem. The MiniBooNE collaboration proposed an effective parameter k 
to increase the effect of Pauli blocking [14 1. CCQE fits were done simultane- 
ously to Ma and k, treated as free parameters. In the recent MiniBooNE's 
paper [9] the best fit for k is within 1 a consistent with k = 1 (means no 
modification of the Pauli blocking). It is important that the 1-parameter fits 
for Ma (with n = 1) do not lead to significantly different results. Also the 
MINOS collaboration proposed an ad hoc modification of the Pauli blocking 



An interesting theoretical idea to explain the Ma value discrepancy comes 
from the sophisticated many-body nuclear model proposed about 10 years 
ago by J. Marteau |15| and developed recently by M. Martini [16]. This is 
the non-relativistic model that includes QE and A production primary in- 
teractions, RPA corrections and local density effects. The most interesting 
feature of the model is the evaluation and inclusion of elementary 2p-2h ex- 
citations. This contribution goes beyond the IA and is absent in free nucleon 
neutrino reaction. There is the strong evidence from electron scattering data 
(even off light nuclei like A He) that in the transverse response function the 
contributions from one and two body mechanisms are similar in strength 
[17j . The 2p-2h contribution is quite large and it is claimed to be able to 
account for the large CCQE cross-section as measured by the MiniBooNE 
collaboration, In the case of neutrino-carbon CCQE process after averag- 
ing over the MiniBooNE beam, nuclear effects are expected to increase the 
cross-section from 7.46 to 9.13 (in the units of 10~ 39 cm 2 ). This includes a 
reduction of the cross-section due to RPA effects and increase thanks to 2p- 
2h contribution. In the case of antineutrino-carbon CCQE reaction the RPA 
and 2p-2h effects cancell each other approximately, and the cross-section is 
virtully unchanged (modification from 2.09 to 2.07). We see that it will 
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be very important to compare predictions of the model with anti-neutrino 
CCQE MiniBooNE data. It will take about 1-2 years before MiniBooNE 
CCQE anti-neutrino data analysis is completed [18J. It is also important 
to compare predictions from the model with MiniBooNE's CCQE double 
differential cross-section data for the distribution of events in muon kinetic 
energy and scattering angle [9] . 2p-2h contribution is parametrized in terms 
of energy and momentum transfer (there are 2 different parametrizations) 
which can be translated into double differential cross-section. The hope is 
that when the extra contribution is superimposed on the IA QE events with 
standard Ma ~ 1-05 GeV the overall distribution mimics the pure QE one 
with M A ~ 1.25 GeV. 



3. Neutral current 7T° production 

Neutral current tt° production (NC7r°) is a dangerous background to u e 
appearance oscillation signal in Cherenkov detectors and during the last 
~ 5 years there were several attempts to measure its cross-section. Since we 
are interested in 7r° leaving nucleus, the experimental data is always given 
in this format with all the FSI effects included. If this definition of NC7r° 
is adopted the relevant events can origin from: NCl7r° primary interaction 
with 7T° not affected by FSI, NCl7r + primary interaction with tt + being 
transformed into 7r° in charge exchange FSI reaction, ... It is clear that a 
comparison of theoretical models with such results is rather difficult. On the 
other hand, such format of the data is very useful to check the performance 
of MC generators of neutrino interactions and tools like GiBUU [19J. An 
additional challenge related with the NC7T production is that it includes a 
COH contribution, which in the case of MiniBooNE's beam neutrino-carbon 
reactions (< E v >~ 1.2 GeV) is estimated to be on the level of ~ 20% [20] . 

Four recent measurements of NC-zr (K2K [21j, MiniBooNE neutrinos, Mini- 
BooNE antineutrinos |22| , SciBooNE |23| ) are complementary. They use 
three different beams (K2K, FermiLab Booster neutrino and anti-neutrino 
beams) and the targets: H 2 (K2K), CH 2 (MiniBooNE) and C 8 H 8 (Sci- 
BooNE). MiniBooNE presented the results in the form of absolutely nor- 
malized cross-section while K2K and SciBooNE reported only the ratio 
da(NClir°)/a(CC). There is an important difference in the measured quan- 
tity: K2K and MiniBooNE present their results as measurements of the final 
states with only one ir° and no other mesons. SciBooNE speaks about states 
with at least one tt° in the final state i.e. a contamination from l7r°l7T , 2-7T 
and > 2ir (with > l-zr ) final states is included (evaluated with NuWro to be 
on the level of 17%). All the experimental groups present the final results 
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in the form of flux averaged distributions of events in the 7r° momentum, 
and in the case of MiniBooNE and SciBooNE in the ir° production angle. 
Additionally, MiniBooNE provides also the data for the cross-section before 
subtraction of the contribution from the neutrino or anti-neutrino beam con- 
taminations. This data represent the measurement that is in the maximal 
possible degree independent on the assumptions of the MC events generator. 

The comparison with the NC7r° data can only be done within MC events 
generators. There are several ingredients of the MC that are tested simul- 
taneously: NC A (or more generally: resonance) production, nonresonant 
contribution to SPP, modification of the A width in the nuclear matter, an- 
gular distribution of the 7r° from a A decay, a COH component, a multipion 
production, FSI effects (the n° absorption rate in the nucleus, cross-sections 
for pion charge exchange reaction in the nucleus, formation zone effects). 
It is an unfortunate situation that the separate ingredients of the MC are 
known with unsatisfactory precision. 

There is only one reported measurement of the NCl7r° cross-section in the 
Gargamelle bubble chamber (the target was composed of C^Hg (90%) and 
CF^Br (10%)) with the ~ 2 GeV beam of neutrinos [23] (we notice that 
the data from the paper [24J contain also a non-negligible contribution from 
the COH reaction what makes the extraction of the resonance NC7r° pro- 
duction even more complicated). The results were published in the form 
of efficiency corrected numbers of events in several exclusive SPP channels 
and also as nuclear effects corrected numbers of events on free nucleon tar- 
gets. The way of treatment of the nuclear effects for the topologies of final 
states is described in [25J . Later reanalysis of the data introduced absolute 
normalization (it was possible because in the original paper one CC SPP 
channel was included in the analysis) and flux averaged cross-sections were 
obtained [26J. An interesting feature of the final results is that the proton's 
NC-7T cross-section is much larger than the neutron's one. Many details of 
the analysis are given in the original paper and it seems possible to reex- 
amine nuclear effects using the better knowledge of pion absorption rates in 
various nuclei. 

The NC7T differential cross-section for pion production angle is very sen- 
sitive to the COH component that contributes to the forward directions 
only. The differential cross section for pion kinetic energy is very sensitive 
to FSI effects and in particular to the absorption rate dependence on pion 
momentum. A useful comparison of performance of FSI models in various 
MCs is done in [27]. Authors of intranuclear cascade codes confront the 
models with the available tt + nucleus scattering data [28J. There is quite a 
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lot of 7r + — 12 C reaction (with the separate absorption and charge exchange 
rates) data in the tt + kinetic energy range 50 — 500 MeV and the typical 
uncertainty is ~ 20% [29]. The data for ir + — ie O absorption is scarse and it 
is not obvious how the reaction cross-sections scales as the size of a nucleus 
increases. Other observables that are useful in benchmarking FSI models 
are nucleon's [30] and pion's [31] [32] transparencies. 



4. Coherent pion production 

In the coherent pion production (COH) the target nucleus remains in the 
ground state. Four possible channels are possible, for CC and NC reactions, 
for neutrinos and for anti- neutrinos, see Eq. @. The NC COH cross-section 
is believed to contribute significantly to the overall NC-n -0 cross-section and 
the process has been a subject of intensive experimental and theoretical 
investigation. There is a clear experimental signal for the COH reaction 
at higher energies and the aim of recent measurements is to fill a gap of 
the knowledge of about ~ 1 GeV COH cross-section. In the case of CC 
reaction K2K [33] and SciBooNE [M| reported no evidence for the COH 
component. In the case of NC reaction MiniBooNE [201 and SciBooNE [35^ 
detected the COH component. The identification of the COH signal is done 
based on the predictions of the MC generatotor of events. In the case of CC 
COH reaction the sample of 1tt + events was searched at low Q 2 , and for the 
NC reaction the distributions of observed it in the forward direction was 
analysed. 

It is interesting to look at the numerical values obtained in the measure- 
ments. For the K2K's analysis of the CC COH reaction, the 90% confidence 
limit for the upper bounds of the COH cross-sections on carbon was esti- 
mated to be 0.6% of the total CC cross-section. Similarly, the SciBooNE's 
results (also for the carbon) are: 0.67% at < E v >~ 1.1 GeV, and 1.36% at 

< E v >~ 2.2 GeV. For the NC reaction, thanks to the information about 
recoil protons, SciBooNE [35] evaluated the ratio of the COH NC7r° pro- 
duction to the total CC cross-section as (1.16 ±0.24)%. SciBooNE reported 
the measurement also in the form of the ratio of CC COH ir + to NC COH 
7T° production, estimated as 0. 14^28 • The value is surprisingly low and is 
a challenge for theoretical models. For the NC reaction MiniBooNE also 
evaluated the COH component (plus possible hydrogen diffractive contri- 
bution about which little is known) in the NC7T production as 19.5% (at 

< E u >~ 1.2 GeV) and then the overall flux averaged NC17T cross-section 
as (4.76 ± 0.05 ± 0.76) • 10- 40 cm 2 /nucleon. Unfortunately, it is difficult to 
translate both measurements into the absolutely normalized value of the NC 
COH cross-section. The reason is that the first from the above two measure- 
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merits strongly depends on the details how is NC7r° production modelled 
in the MC used by the MiniBooNE collaboration i.e. NUANCE [11] • In 
NUANCE, RES, COH and BGR (background) NCvr reactions are defined 
according to the primary interaction. A peculiar feature of NUANCE is that 
COH produced pions are subject to FSI, in a similar way as pions produced 
in the RES primary interaction. In the final step of the analysis the fit is 
done for the composition of the sample of NC7r° events in terms of three 
components, and the COH fraction is defined as xqoh I '{xcoh + xres) [57] - 

NUANCE contains an independent diffractive (DIFFR) component of SPP 
coming from reactions on free nucleons [36]. The experimental signal of 
pions produced by the DIFFR mechanism is similar to the COH one, and 
in the process of selection of COH events they are likely to be put to the 
same category. In |36j the evidence is shown for the presence of DIFFR 
component in the high energy neutrino reactions, in the region of invari- 
ant hadronic mass W > 2 GeV. For W < 2 GeV the DIFFR component 
is expected to contribute to the nonresonant background and in the MC 
implementation it is important to avoid a danger of double counting. Mini- 
BooNE estimated the DIFFR contribution as 16% of the overall 'coherent' 
cross-section but its impact on the reported value of the COH fraction is 
not very important. 

The evaluation of the COH contribution to the cross-section is always made 
by confronting the measurements with the predictions of MC events gener- 
ators. Before the last Nulnt09 workshop S. Boyd and S. Dytman initiated 
the theory and MC comparison project [38]. In the case of COH reaction 
ten different computations were compared, four from MC generators (GE- 
NIE [40], NEUT [41], NUANCE, NUWRO 021) and six from theoretical 
models. All the MC generators implemented the same theoretical model of 
Rein&Sehgal [43] . It was surprising that the predictions from the MCs differ 
by more than 100% . For the CC COH reaction on carbon at E v = 1.1 GeV 
the predicted cross-sections (per nucleon in the units of 10 _40 cm 2 ) were: 
2.33 (NEUT), 1.88 (NUANCE), 1.31 (NUWRO) and 0.85 (GENIE). Pre- 
dictions from theoretical models were usually about 0.5, slightly larger were 
results from the computations of models of Alvarez-Ruso et al [44] (0.8) and 
Martini et al [16] (0.66). 

It is sometimes believed that it is sufficient to correct the implementation of 
the model [43] by a correction factor (e.g. 2/3 as it was done by MiniBooNE 
in the case of NUANCE) . But the comparison of differential cross-sections in 
pion kinetic energy revealed even more dramatic differences between MCs 
|38j . Theoretical computations predict smooth distributions with a max- 
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imum at T n o ~ 150 MeV (the actual position varies from 110 MeV to 
180 MeV). The shape is correctly reproduced by NUANCE, while GENIE, 
NEUT and NuWro predict a lot of extra structure: a local minimum at 
T^o ~ 170 MeV (presumably from pion absorption, GENIE, NuWro) or at 
T n o ~ 300 MeV (NEUT). The large differences between MCs and theoreti- 
cal models are clearly seen also in double differential cross-section at fixed 
pion production angle. 

Higher neutrino energy {E u >~ 2 GeV) COH production (including re- 
cent NOMAD measurement) was successfully explained with the help of the 
PCAC based model [13]. Adler's theorem relates (Jcoh{v+X — > v+X+tt ) 
at Q 2 ->• with cr(vr° + X ->• vr° + X). Recently the model for the CC re- 
action, has been refined |45j to include charged lepton mass effects. The 
new model predicts the <tcoh(^ + ) I '^COh{^) ratio at E v = 1 GeV to be 
1.45 rather than 2. Another very important modification was to use the 
available data for da(-K + 12 C — > it + 12 C/dt in the region of pion kine- 
matical energy 100 MeV< T^ < 900 MeV. As a result the predicted COH 
cross-section becomes smaller by a factor of 2-3. The other PCAC based 
approach is proposed in |46j . The microscopic models for the COH reaction 
|47j |44| |16| assume A dominance and are expected to be more reliable 
at low neutrino energy. Within microscopic models there are still various 
approaches e.g due to differences in the treatment of a nonresonant back- 
ground. The absolute normalization of the predicted cross-section depends 
on the adopted value of the TV — > A form factor C§ (0) . 

5. Other measurements 

For completness I list and comment some other recent neutrino cross- 
section measurements and phenomenological investigations. 



MiniBooNE collaboration measured the ratio a(CClir + )/a(CCQE) 
which is an important observable, independent on the overall neutrino flux 
normalization uncertainty. In the MC independent version the reported ob- 
servable is the ratio of SPP-like to QE-like cross-sections as a function of 
neutrino energy. The uncertainty of this measurement is still larger than 
10%, and the reason is the error in the recontructed neutrino energy. Us- 
ing information contained in MC model it was also possible to evaluate the 
ratio of SPP to QE defined as primary interactions on an abstract isoscalar 
nuclear target. From the MiniBooNE data it follows that in the Monte 
Carlo events generators an increase of the axial mass Ma, which controls 
the value of the CCQE cross-section, must be accompanied by a change of 
parameters that determine the value of CCl7r + production. 



12 epifania'acta'15'maja printed on May 20, 2010 



MiniBooNE reported preliminary results for the CCl7r + -like cross-section 
|49j . Thanks to the possibility to detect final n + , it was possible to recon- 
struct neutrino energy and invariant hadronic mass W, based on muon and 
pion observables only. W distribution shows a clear A peak. Both the total 
and differential cross-sections disagree with the MC (NUANCE) predictions 
and the data can be very useful for the future improvements in the numer- 
ical codes. 

Almost all the MC events generators rely on the Rein-Sehgal model for pion 
resonance (RES) production |50j. The model includes contributions from 
18 resonances and covers the region W < 2 GeV. In the original model, 
the charged lepton is assumed to be massless. The RS model should be 
modified in order to be more reliable in the sub-GeV neutrino experiments. 
Lepton mass corrections can be introduced following the prescriptions pro- 
posed in [51] . There are also arguments in favour of modifying form- factors 
used in the model [52]. It is known that the vector part does not reproduce 
the electron scattering data and consequently also the axial part should be 
fixed by making a fit to the deuteron and/or MiniBooNE pion production 
data. In the paper [53] a comparison study of some theoretical models [52] 
with the MiniBooNE data was done. 



MiniBooNE reported also the results with NC elastic reaction cross-section 
[56]. The measurement was possible because the MiniBooNE detector, 
which is basically a Cherenkov detector, can observe also scintillation light 
from low momentum nucleons. An attempt is done to measure the strange 
quark component of nucleon spin using the proton enriched sample of events. 

6. Monte Carlo generators 

The performance of MC events generators cannot be better than the 
precision of neutrino cross-section measurements and this is on the level of 
15%-20%. This explains the differences in the performance of MCs |27j . 
Most important MCs (GENIE [3D], NEUT [UJ, NUANCE [llj , FLUKA 
[57| ) are principally aimed to be tools in the data analysis and their devel- 
opments are guided by the needs of particular experiments [58]. Each of 
them has some priorities determined by the physical program and by the 
detection techniques. From the experimental point of view it is useful that 
the performance of a MC can be fine-tuned to a particular experimental 
situation and that one can rescale contributions from some exclusive chan- 
nels in order to get an agreement with what is observed. From time to time 
there are discussions (e.g during the first NuIntOl, and also during the last 
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Nulnt09 [59J) about the construction of a universal Monte Carlo generator 
of neutrino interactions. Few years ago a project aimed in this direction 
was launched under the name of GENIE {Generates Events for Neutrino 
Interaction Experiments) |40j. It is written in C++ and its architecture 
allows to absorb various pieces of codes from other MCs. In the core of 
GENIE there is another MC generator NEUGEN [60], and there is still an 
intensive work on various elements of the code. 

It seems however that in the perspective of 5 years several MCs will be 
developed in parallel, what is probably useful for their quality: it will be 
possible to compare and cross check the performance of proposed improve- 
ments. In the past many MC comparison tests were done and presented 
during Nulnt workshops [2J. Currently, in the T2K collaboration, several 
useful GENIE-NEUT comparison studies are performed. It is also impor- 
tant that there are MCs or MC-like tools developed by theorists (GIBUU 
|19j . NuWro [12]) which allow for more flexible manipulations with various 
ingredients of the codes. MCs are large and long-term projects and it is 
often difficult to verify quickly all the consequences of even minor modifica- 
tions in order to avoid risks of spoiling the self-consistency of the code. 

There are very few obvious modifications which definitely should lead to 
an improvement of the performance of MCc. One of them is, as it was 
argued earlier, an implementation of a better model for the COH pion pro- 
duction. Another one seem to be the replacement of the Fermi gas (FG) 
model by a more reliable formalism in the description of a nucleus target. 

Fermi gas model is determined by two parameters: Fermi momentum and 
binding energy and it defines the probability distribution of finding inside 
nucleus a nucleon with a given value of momentum (quadratic distribution) 
and binding energy (a constant value). From the electron scattering exper- 
iments it is known that in the interesting kinematical region the FG model 
allows for reasonable agreement with the data. The advantage of the model 
lies in the simplicity. Also its MC implementation is straghtforward. Is is 
however known from the more detail analysis of the electron scattering data 
that the FG model suffers from many limitations and is unable to reproduce 
separately longitudinal and transverse nuclear response functions. 

There is a variety of approaches to describe nuclear effects in electron scat- 
tering which were later applied to neutrino interactions. Short descriptions 
of many of them is given in [38] [39] together with a list of basic refer- 
ences. Monte Carlo implementation of most of the theoretical models is 
usually complicated, if not impossible. MC requires a propagation of all 
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the degrees of freedom (particles taking part in interactions) and it is not 
sufficient to know the inclusive neutrino interaction cross-section. There is 
however a model that can replace FG and which allows for much better 
agreement with the data. It is called spectral function (SF) [61] but strictly 
speaking the name hole's spectral function should be used. The use of SF 
has been advocated by Omar Benhar during several Nulnt workshops. SF 
is defined as a joint probability distribution to find inside nucleus a nucleon 
with a given momentum and binding energy. SF arises naturally as one 
calculates the neutrino quasielastic cross-section in the plane wave impulse 
approximation (PWIA) [62] i.e. assuming that the nucleon arising in the 
final state after primary interaction leaves the nucleus with no FSI effects. 
Together with a model for FSI [63j (or for the particle spectral function) 
the SF model leads to a very good agreement with the electron-nucleus 
cross-section data for momentum transfers larger than ~ 350 GeV [64] . 
The available models of SF combine information from the mean field theory 
(shell model) and a contribution from short range correlations (SRC). The 
shell model structure is clearly seen as the peaks in the probability distribu- 
tion corresponding to shell model orbitals. SRC part contributes to a large 
nucleon momentum tail in the probability distributions. SF for few impor- 
tant nuclei (carbon, oxygen, iron) were constructed by Omar Benhar and 
his collaborators. There also exist approximate models of SF for medium 
size nuclei like calcium and argon, which were shown to lead to a good agree- 
ment with the electron scattering data [55]. It is estimated that correlated 
proton-neutron pairs are ~ 18 times more likely than proton-proton ones, 
and that ~ 25% of nucleons inside medium size nuclei have momenta larger 
than 300 MeV [66]. Neutrino QE interactions may provide an opportunity 
to observe correlated pairs because the signal (up to distortions caused by 
FSI effects) would be a pair of high momentum protons possible to observe 
e.g. in liquid argon detectors. 

The implementation of SF formalism into MC events generators is not triv- 
ial. SF introduces extra integrations which, if not performed in the numer- 
ically efficient way, can make the process of generating events very time 
consuming. For that reason a few years ago an effective approach was pro- 
posed [67] . The genuine SF defines the momentum distribution and also the 
momentum dependent binding energy, which is defined so that it replaces 
the constant binding energy of the FG model. It was shown that the effective 
approach provides a good approximation of the SF [67j . Its implementa- 
tion is easy unless the MC code (that is the case for NEUT) uses analytical 
formulas obtained by averaging elements of the hadronic tensor over target 
nucleon momentum within the FG model. Recently, the SF formalism was 
implemented in the MC generator NuWro [68]. It has been already known 
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that SF does not introduce much change, with respect to FG, as far as the 
shape of the QE differential cross-section in Q 2 is concerned [69]. In [68] 
it is shown that the fit to the value of the axial mass Ma based on the 
MiniBooNE's double differential cross-section in muon kinetic energy and 
scattering angle done within the SF model leads to a very similar value as 
the FG model (Ma ~ 1.4 GeV). Similar results were obtained also by the 
MiniBooNE collaboration [70]. In the case of FG model the goodness of fit 
is better because the use of SF reduces the predicted value of overall CCQE 
cross-section by about 20%. The elimination of bins with large contribu- 
tion of events with small values of momentum transfer (more than 50% of 
events with |g| < 350 MeV) makes the fitted value of Ma smaller by about 
100 MeV. The modification goes in the desired direction but there remains a 
large gap with respect to the old deuterium bubble chamber measurements. 
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